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Measurements of gamma rays with Nal(Tl) crystals yield pulse height distributions 
related to the true energy spectrum by an integral equation whose kernel (response function) 
is the probability that an incident photon of energy E will give rise to a pulse of size E' . 
The response function has been calculated by the Monte Carlo method for photons with 
energies from 0.279 to 4.45 Mev. Results are presented for cylindrical crystals ranging in 
size from 0.25 (radius) by 0.5 inch (length) to 2.5 by 9 inches. They are based on the 
analysis of 50,000 photon histories sampled with the use of the Standards Electronic Auto- 
matic Computer (SEAC) . Analytical corrections were made for the escape of annihilation 
radiation and bremsstrahlung from the crystal. 



1. Introduction 

The response function of a Nal(TI) crystal spec- 
trometer relates the observed pulse height distribu- 
tion to the true gamma-ray spectrum. The primary 
means of determining this function is, and pre- 
sumably will remain, experimental. But a good 
theoretical understanding will contribute toward 
improving scintillation spectrometry. 

While the theory of the diffusion of gamma radia- 
tion in an infinite medium is now well developed, the 
calculation of the response function is a boundary 
problem whose exact solution is beyond the scope of 
theory in its present state. 

The subject of this article is a Monte Carlo calcula- 
tion of the response function, i. e., an experiment 
carried out on paper. While this falls short of an 
analytical theory, it is relatively simple to do, and 
has important advantages compared to a physical 
experiment. One is not hampered by the limited 
availability of monoenergetic sources nor by acci- 
dental disturbing effects such as background radia- 
tion. 

Exploratory Monte Carlo calculations of the 
response function have already been made by Camp- 
bell and Boyle l for energies E^6 Mev, and by Foote 
and Koch 2 for E^4A5 Mev. An approximate 
analytical calculation for E^l.3 Mev was carried 
out by Maeder, Miiller, and Wintersteiger. 3 These 
authors computed the effect of the first interaction 
of the incident photon in the crystal exactly and made 
an elaborate estimate of the effect of multiple inter- 
actions. This approximation limits the validity of 
their considerations to small crystals not much larger 
than one mean free path of the incident radiation. 4 
They confined their attention mostly to cylindrical 
crystals whose radius equals their length. In view 

"This work was supported by the United States Atomic Energy Commission. 

i J. G. Campbell and A. J. F. Boyle, Australian J. Phys. 6, 171 (1953). 

2 R. S. Foote and H. W. Koch, Rev. Sci. Instr. 25, 746 (1954) . 

3D. Maeder, R. Muller, and V. Wintersteiger, Helv. Phys. Acta 27, 3 (1954). 

4 Maeder et al. also attempted a calculation for an "infinitely large" crystal, 
i. e ., for a semi-infinite scattering medium, by considering the first two interactions 
exactly, and estimating the effect of subsequent interactions. As we shall show 
in section 4.1, this procedure was inaccurate (except at low energies) so that it 
does not provide a good basis for extrapolating their results from small to large 
crystals. 



of the complexity of their "orders of scattering'" 
calculation and of the approximations whose con- 
sequences cannot be easily surveyed, a direct re- 
calculation of their results seemed of value. More- 
over, one would also like to know more about the 
influence of the crystal shape on the response 
function. 

Monte Carlo calculations have been made of the 
response functions of cylindrical crystals of diverse 
shapes ranging in size from small to very large for 
radiation incident with energies of 0.279, 6.661, 1.17, 
1.33, 2.62, and 4.45 Mev. 

The absorption and Compton scattering of photons 
in the crystals were calculated by random sampling, 
the computations being carried out on the SEAC. 
Approximate analytical corrections were applied to 
the Monte Carlo results to account for the escape of 
secondary bremsstrahlung and annihilation radiation 
from the crystals. 5 These corrections were quite 
small at 1.17 and 1.33 Mev, but appreciable at 2.62 
and 4.45 Mev. Because of the approximate nature 
of the corrections, the most important and accurate 
parts of this work are those pertaining to the four 
lowest energies where the corrections were negligible. 
Most of the results are for collimated sources, while 
some sample calculations pertain to broad beams, 
point-iso tropic sources, and off-axis collimated sources. 

The main results of this investigation are: (1) The 
confirmation and expansion of the results of Maeder 
et al. for small crystals; (2) the extension to large 
crystals in which multiple interactions are impor- 
tant and to higher energies; (3) last but not least, 
the creation of a SEAC code by means of which re- 
sponse functions for source and crystal geometries 
peculiar to a given experimental situation can be 
readily computed with a modicum of effort. 



5 One of the reasons for this approximate procedure was the temporary non- 
availability of SEAC during the latter part of this investigation. It would 
have been better to carry through the entire calculation on the automatic com- 
puter. This is particularly easy for quanta resulting from the annihilation of 
positrons produced in a first interaction of the primary photon. The source 
distributions of these photons within the crystal can easily be calculated, and 
their subsequent history traced by means of the existing SEAC code. We hope 
to change the SEAC code in the future so that it can handle the escape of all 
secondary radiation, and to extend the calculations above 4.45 Mev. In the 
meantime, there was some merit in reporting the present results for 2.62 and 4.45 
Mev, since calculations in this energy range have not been reported previously. 
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2. Method of Calculation 
2.1 Formulation of the Problem 
The response function can be expressed in the form 



R(E,E")= rL(E,E')g(E',E")dE' 



(1) 



where L(E,E')dE' is the probability that a photon 
of initial energy E will deliver to the crystal — as the 
result of one or more interactions — an amount of en- 
ergy between E 1 and E' + dE') and g(E',E")dE" is 
the probability that upon delivery of energy E f to 
the crystal, the light output as amplified by the 
attached photomultiplier will indicate an apparent 
energy (pulse height) between E n and E f,J rdE'\ 
It is an experimental fact that 

g(E' ,E")=(2ttcxET 1/2 exv{-(E'-E") 2 /2aE'}. (2) 

The value of the parameter a depends on the physical 
characteristics of the crystal-photomultiplier system. 
We shall consider this parameter as determined ex- 
perimentally, and concentrate on the calculation of 
the loss function L(E,E'). 

To clarify the physical picture, it is useful to break 
the loss function up into a number of component 
parts, writing 

L(E,E')=Y(E){K(E,E')+p{E)5(E-E')}. (3) 

Y(E) is the efficiency, i. e., the probability that a 
photon incident with energy E will have at least one 
interaction (scattering or absorption) in the crystal. 
For a collimated source 



Y(E)=l-exp [-n(E)L] 



(4) 



where L is the length of the crystal, and ix(E) the 
narrow beam linear absorption coefficient. 

K(E,E')dE' is the probability that an interacting 
photon will emerge from the crystal after one or 
more Compton scatterings, leaving behind an 
amount of energy between E' and E'-\-dE'. The 
function p(E) represents the probability that an 
interacting photon will be absorbed (possibly after 
a number of Compton scatterings), while the delta 
function 8(E—E') indicates that the entire energy 
E is delivered to the crystal in a photon history 
terminating in absorption. We note the normaliza- 
tion 



s: 



dE'K(E,E')+p(E)=l. 



(5) 



For a wide range of conditions the shape of the loss 
function depends largely on the ratio of the numbers 
of absorbed and emergent photons, and only rather 
insensitively on K(E,E'). In a rough approxima- 



tion, the shape can therefore be characterized by 
p(E), which — following Maeder — we shall call the 
photofraction. In experimental terms, 



V(E)- 



area under the "photopeak" of the pulse 

height distribution 
area under the entire pulse height distribution 

(6) 



Knowledge of the photofraction is useful in other 
connections. Thus if one wants to use a crystal for 
counting photons in the presence of a high back- 
ground, it may be desirable to count only the large 
pulses resulting from complete photon absorption. 
The required "photoefficiency" of the crystal is 
Y(E) p{E). Counting of photons in the "photo- 
peak" only is particularly useful for disentangling 
complex spectra, containing several lines. 6 

2.2. Random Sampling 

For the calculation of the loss function, a direct 
stochastic analog method paralleling the physical 
processes in all respects, has been chosen. Scintilla- 
tion crystals are very efficient detectors; hence the 
efficiency of an analog Monte Carlo calculation is 
also high, and the computational cost of statistical 
and analytical refinements would be out of proportion 
to the possible gain in efficiency. 

The calculation proceeds along the following lines. 
We start a photon at a specified position on the 
crystal surface with specified energy and direction. 
The position of the first interaction is determined 
as if the photon moved in an infinite scattering 
medium ; the distance traveled to the point of inter- 
action is distributed exponentially, the mean of the 
distribution being the inverse of the total linear 
narrow beam absorption coefficient. If the inter- 
action occurred outside the region occupied by the 
crystal, the photon is recorded as having had no 
interaction, and another photon history is begun. 
If the interaction is inside, we determine by random 
sampling whether it is an absorption or Compton 
scattering, the relative probabilities for these contin- 
gencies being proportional to the respective cross 
sections. In case of an absorption, the photon 
history is terminated, and another begun. In case 
of a scattering, the new energy and direction of the 
photon after the collision are sampled from the 
Klein-Nishina distribution, whereupon the location 
of the next interaction is sampled in the same manner 
as before. If this interaction occurs outside the 
crystal, the photon is considered to have escaped, 
and the energy left behind in the crystal is recorded. 
If the interaction was inside, we proceed as before, 
deciding whether it was an absorption or scattering, 
etc. The procedure is repeated until the photon 
has either left the crystal, or has reached an energy 

W. E. Kreger, Phys. Rev. 98, 1554 (1954); see also footnote 3. 
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lower than 50 kev, in which case the photon is 
considered as absorbed. The bias introduced by 
this cutoff is negligible because of the steep rise of the 
photoelectric absorption cross section for Nal at 
energies below 50 kev. 

The detailed equations for carrying out the 
sampling procedure have been presented elsewhere, 
together with an account of the numerical manipula- 
tion required to adapt the various cross sections to 
computation on a high-speed computer. 7 

The required absorption coefficients were obtained 
from the tabulation of G. White. 8 

2.3. Machine Computation 

In the SEAC program, the calculation was set up 
for the simultaneous consideration of nine crystals 
of different sizes. For the sake of. convenience, these 
sizes were chosen so that the crystals formed a nested 
sequence, each crystal containing completely the next 
smaller one. When a photon in the random sampling 
experimenl escaped from a crystal, the history was 
continued for all crystals of larger size. This pro- 
cedure resulted in great economy of computing time. 9 
Beyond that, it increases the statistical efficiency of 
comparisons of the response functions of crystals of 
different shapes. 

The output of a SEAC calculation consists of the 
following information for each of nine crystals: 

a. The efficiency Y(E). This information is in a 
sense redundant since it can also be easily obtained 
by a direct analytical calculation (eq (4)). But a 
comparison between analytical and Monte Carlo 
efficiencies provides a useful check that the computer 
operation was error-free. 

b. The spectrum (in histogram form.) of the energy 
packets supplied to the crystal by photons which 
eventually escape: 



n+J 

UE)= f 32 i 
J "e 



K(E,E')dE' (11=0,1, . . • 31). 



c. The photofraction p{E). 

It was convenient to generate and process photon 
histories in groups of 1,C00. For photons with an 
initial energy of 660 kev, this took about 35 to 40 min 
of SEAC time; for higher energies the running time 
was a little longer. 

2.4. Escape of Secondary Radiation 

In the Monte Carlo calculations it was assumed 
that the entire energy released in the successive 
scatterings or final absorption of the incident photons 

' M. Berger, J. Research NBS 56, 343 (1956) RP2640. 

8 G. White (private communication). 

• There is another obvious way of increasing the efficiency of the computation, 
namely, by treating the first interaction analytically, and sampling only subse- 
quent interactions. This procedure avoids photon histories which result in no 
interaction at all, and thus contribute only information on the efficiency (which 



remained in the crystal. Actually this is the case 
only at low energies. At energies above the threshold 
for pair production, annihilation quanta may escape 
from the crystal; additional energy can be lost in the 
form of bremsstrahlung. In the energy region of 
interest (up to 4.45 Mev) the effect of bremsstrahlung 
was found to be minor compared to that of annihila- 
tion radiation. Not only is the latter greater 
quantitatively, but it also affects the shape of the 
response function in a more distinct fashion, giving 
rise to two additional peaks corresponding to the 
complete escape of 1 or 2 annihilation quanta 
subsequent to the absorption of a primary photon. 
The appropriate corrections were made according 
to the approximate procedure outlined in the 
appendix. The photofraction p(E) was corrected 
both for annihilation radiation and bremsstrahlung. 
In order to avoid excessive computation, the detailed 
shape of the response function was corrected only for 
annihilation radiation, and an adjustment for the 4 
effects of bremsstrahlung was made by a subsequent 
renormalization of the function K(E,E').- 

2.5. Gaussian Broadening of the Response Function 

It follows from eq (1), (2), (3), and (4) that the 
response function ^j 



R{E,E")= 



r,E) 

-\2ttolE 



$° K{E > E,) w^ { - {E ' 

~EJI2aE'}dE' 



Y(E)p(E ) 
■y/2iraE 



exp {-{E"-Efj2aE} 



y (7) 



*H(E,E")+G(E,E") 

(tail) (photopeak) 



For numerical applications the parameter a was 
chosen to be 0.0014G E(Mev), corresponding to a 

width at half-maximum W=0.09^jE. This is a 
typical value in close agreement with experimental 
data of Foote and Koch (see footnote 2). 

In the evaluation of H(E,E") one must take into 
account that K(E,E') is available only in histogram 
form 

71+1 



*-/. 



K{E,E')dE', 



71 = 0,1, 



31. 



— E 
32 



in any case can be determined analytically), but not the shape of the loss function. 
This procedure leads to complications when one tries to consider several crystals 
simultaneously, and it was decided that— on balance— the possible gain in 
efficiency was not sufficient to compensate for the increased amount of manipula- 
ton and computing time. 
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It is assumed that the k n 's have been corrected if 
necessary, for escape of secondary radiation. We 
consider them to form a row matrix k(E) = (kJci . . . 
&31). We similarly replace H(E,E") by a row matrix. 



h(E)=(hohi . . . h 3 i), where h 



• (e) =j. 



32~ J 



H(E,E')dE'. 



32 



The convolution of K(E,F/) with a Gaussian is 
replaced by a matrix multiplication: 



h=TL 



(8) 



where the matrix element T nm (a kind of transition 
probability) is given by 



T nm = dt 



j: 



1 



J2wa 



n+t 
32 



ro+1 

dE' exp /-(#' 

32 



(9) 



The exact location of a pulse within an interval of 
size E/32 is not known from the Monte Carlo cal- 
culation. We assume a uniform distribution, which 
is the reason for the average with respect to t in (9). 



3. Results 

3.1. Efficiency and Photofraction for Collimated 
Radiation 10 

Crystals of such sizes as are either readily avail- 
able commercially, or are in use in this laboratory, 
were selected for computation. In table 1, the crystal 
efficiencies are listed for collimated radiation of 
various energies incident along the axis of the cylin- 
drical crystals, computed according to eq (4). The 
figures in parentheses indicate the values that were 
obtained by the Monte Carlo calculations. The 
agreement is excellent. 

In table 2, the photofraction p(E) is given. This 
is the corrected value taking into account the escape 
of secondary radiation. The amount Ap which had 
to be subtracted from the raw Monte Carlo result 
p in order to make the correction is indicated in 
parentheses. The indicated errors are standard 
deviations computed as follows. The standard 
deviation of p is 

~L YN 



0o= 



1 



(10) 



where N is the total number of photon histories con- 
tributing to the determination of p . (N is also 
shown in tables 1 and 2.) The correction Ap was 
assumed to have a standard deviation (l/10)Ap. 
This is a nominal figure representing an "educated 
guess. " The total standard deviation is 



tr=[oS+(l/10Ap) : 



11/2 



(11) 



i° Collimated radiation means a narrow pencil of radiation incident along the 
crystal axis. 



Table 1. Crystal efficiency 

The values listed without parentheses were calculated according to eq (4); those in parentheses are the corresponding Monte Carlo results obtained through the 
analysis of the indicated number of photon histories. 



Crystal Xo. 
Radius (in.) 
Length (in.) 



Energy E 
(Mev) 

0.279 
.661 
1.17 
1.33 
2.62 
4.45 



1 
2.5 
9.0 



2 


3 


4 


5 


6 


7 


8 


2.5 


2.5 


0. 875 


0.75 


0.75 


0.75 


0.5 


8.0 


4.0 


2.0 


2.0 


1.5 


1.0 


1.0 



0.25 

0.5 



Efficiency Y(E) 



000 
000) 



(. 986) 



958 
(. 964) 



f .950 
I (-947) 



1.000 
(1. 000) 


0.998 
(. 999) 


0.958 
(.957) 


0.958 
(.957) 


0.908 
(. 903) 


0.796 

(. 798) 


0.796 
(. 798) 


0.996 
(. 996) 


.936 
(. 938) 


.747 
(. 748) 


.747 
(• 748) 


.643 
(. 644) 


.497 
(.497) 


.497 
(.497) 


.981 
(. 979) 


.861 
(. 859) 


.627 
(. 635) 


.627 
(. 635) 


.522 
(. 529) 


.389 
(. 395) 


.389 
(. 395) 


.975 

(• 977) 


.841 
(. 833) 


.602 
(. 599) 


.602 
(. 599) 


.499 
(. 495) 


.369 
(. 363) 


.369 
(. 363) 


.940 
(. 946) 


.755 
(. 760) 


.504 
(. 507) 


.504 
(.507) 


.410 
(. 409) 


.296 
(. 292) 


.296 
(. 292) 


.930 
(.927) 


.735 
(• 741) 


.486 
(. 484) 


.486 
(.484) 


.393 

(. 398) 


.283 
(. 282) 


.283 
(. 282) 



0. 548 
(.546) 

.290 
(.295) 

.218 
(. 223) 

.206 
(• 191) 

.161 

(• 164) 

.153 
(. 153) 



No. of 
histories 

N 



6,000 
10, 000 
5,000 
5,000 
6,000 
5,000 
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Table 2. Photof taction p{E) 

The values listed without parentheses pertain to photofractions corrected for the escape of annihilation radiation and bremsstrahlung. The corresponding figures 
in parentheses indicate the correction A/> thai had to he subtracted from the raw Monte Carlo values Po(E). The indicated standard devial ions were computed ac- 
cording to eq (11). 



Crystal No. 
Radius (in.) 
Length (in.) 



Energy E 

(Mev) 



2.62 



2.5 
9.0 



2 


3 


4 


5 


6 


7 


8 


2.5 


2.5 


0.875 


0.75 


0.75 


0.75 


0.5 


8.0 


4.0 


2.0 


2.0 


1.5 


1.0 


1.0 



Photofraction p{E) 



' 0.973 
±. 002 

' .887 
±. 003 

: .778 
±. 006 
[ (.002) 

767 
±. 006 
(. 003) 

. 653 

±. 008 
(.051) 

.621 

±.015 

(. 139) 



9 

0.25 
0.5 



0.973 
±.002 


0.971 
±.002 


0.915 
±.004 


0.900 
±.004 


0.882 
±.004 


0.843 
±.005 


0.814 
±.006 


.884 
±. 003 


.821 
±.004 


.542 
±.005 


.508 
±.006 


.481 
±.006 


.442 

±.007 


.377 
±.007 


.775 

±. 006 

(. 002) 


.682 

±.007 

(. 002) 


.368 

±.008 

(. 003) 


.343 

±.008 
(. 003) 


.315 

±. 009 
(. 005) 


.282 

±.010 

(. 005) 


.235 

±.010 

(. 005) 


.758 

±. 006 

(. 004) 


. 667 

±.007 

(.005) 


.351 

±.009 

(. 009) 


.325 

±.009 

(. 009) 


.303 

±.009 

(. 010) 


.272 

±.010 

(.013) 


. 224 

±. 010 

(.013) 


. 643 

±. 009 

(. 056) 


. .531 

±.010 

(.069) 


.243 

±.014 
(• H7) 


.226 

±. 015 

(. 122) 


.201 

±. 016 

(• 128) 


.175 

±.017 

(. 135) 


. 144 

±.017 

(.131) 


. 608 

±.015 

(.141) 


.491 

±.020 

(• 184) 


. 169 

±. 033 

(.319) 


.157 

±.033 

(.319) 


.132 

±. 035 
(. 332) 


. 119 

±. 036 

(.331) 


.0826 

±.037 
(.317) 



0.711 

±.008 



.140 

±.011 
(. 003) 

.120 

±.011 

(.011) 

.0950 

±.018 
(. 130) 

.0890 

±.039 
(. 349) 



No. of 

histories 

N 



6, 000 



10, 000 



5. 000 



3.2. An Empirical Formula for the Photofraction 

It has already been noted by Maeder et al., that 
for small crystals at an energy of 0.511 Mev the 
photofraction is a smooth function of ^LR where L 
is the crystal length and R the radius. The results 
of table 2 show that this also holds true for large 
crystals and in the entire energy range from 0.279 
to 4.45 Mev. 

Moreover, the dependence of the photofraction on 
the crystal size can be accurately represented by a 
simple empirical formula, under conditions where the 
correction due to the escape of secondarv radiation 
is absent or small (i. e., at 0.279, 0.661, and 1.17 
Mev for all the crystals, and at 1.33 for all but the 
smallest crystal considered.) 11 

The formula is: 



--A{E)-B(E)e- c ^^ LR (12) 

where L and R are the length and radius of the 
crystal. With the parameters A, B, and C listed 
in table 3, formula (12) represents the values of 
p(E) given in table 2 with an accuracy that is 1 to 2 
percent in most cases, occasionally 3 to 4 percent, 
and in the worst cases 5 to 7 percent. These devia- 
tions are presumably due to the fact that p(E) de- 
pends on L and R in a more complicated way than 
through t/LR. 

ii iAt higher energies a simple expression of the type (12) does not appear to 
work well. In any case the high-energy results are not precise enough to make 
fitting by a formula worthwhile. 



Table 3. Parameters of the empirical formula, eq (12), for 
the photofraction 



E 


HE) 


B(E) 


C(E) 


Mev 






cm~ { 


0. 279 


0. 973 


0. 514 


0. 658 


. 661 


. 940 


. 841 


. 219 


1. 17 


. 980 


. 934 


. 124 


1. 33 


.985 


. 942 


. 120 



Formula (12) also predicts the correct number 
albedo for a semi-infinite medium. When L-^co 
and R->co 

p(*>,E)=A(E)=l— number albedo. (13) 

When we extrapolate in the other direction 
(i-^0 and R-^0) 



2)(0,£) = A(E)-B(^) = mph/m 



(14) 



where ju ph /m is the ratio of the photoelectric to the 
total attenuation coefficient. Thus, of the three 
parameters A, B, and C, two are fixed by (13) and 
(14), and only C, a quantity with the dimensions 
of an attenuation coefficient, is adjustable. This 
makes it rather remarkable that a simple formula 
like (12) works. 
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3.3. Shape of the Response Function for Collimated 
Radiation 

The response functions for all the crystals listed 
in tables 1 and 2 (except No. 2) are shown in figures 
l,a and l,b for incident energies of 0.279, 0.661, 
1.17, 1.33, 2.62, and 4.45 Mev. For the sake of 
clarity of presentation we have plotted only the 
response function exclusive of the photopeak, i. e., 



the function H(E,E') as defined in eq (7). The 
photopeaks are indicated by arrows. 

The Gaussian broadening has been calculated 
according to the prescription given in section 2.5. 
It should be mentioned that the shape of the response 
function is rather insensitive to the magnitude of a. 
The results are presented in histogram form, the 
range between the energy of incidence E, and zero 
energy being covered by 32 intervals of equal size. 
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Figure 1. 



Response functions of Nal(Tl) scintillation crystals for collimated radiation incident along the crystal axis. 

a, Crystals No. 9, 8, 7, and 6. 



Only the "Compton tails" of the spectra are shown; the photopeaks are merely indicated by arrows. All curves are normalized so as to correspond to one incident 
photon the area under them being Y(E) [l-p(E)]. The curves for different energies have been shifted apart through multiplication by indicated scale factors 10a. 
The shaded areas represent the pulses which have been shifted from the photopeak because of the escape of annihilation radiation. 
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The histograms are normalized absolutely so that 
they correspond to one photon incident on a 
crystal. The area under each histogram is thus 
Y(E)[l—p(E)]. 12 In order to supply the missing 
photopeak, one need merely add a Gaussian with 
the desired standard deviation (cf. eq (7)) and 
multiplied by a weight factor Y(E)p(E); Y (E) 
and p(E) can be found in tables 1 and 2. 

The shaded parts of the histograms indicate the 



12 To separate the histograms for different energies, the ordinates were multi- 
plied by indicated scale factors 10a. 



contribution to the scintillation spectrum due to the 
escape of secondary radiation. In other words, they 
represent pulses which would have ended up in the 
photopeak if the absorptions had not been spurious. 
It can be seen that at sufficiently high source energies 
the response functions have two peaks at energies 
mc 2 and 2 mc 2 below the incident energy E, which 
are associated with the escape of 1 or 2 annihilation 
quanta. 

It can be seen that the shapes of the response 
functions for various crystal sizes, and even for 
different energies, are on the whole rather similar. 




8 16 24 32 8 16 24 32 8 16 24 32 8 16 24 32 



»£ 

Figure 1. Response functions of Nal(Tl) scintillation crystals for collimated radiation incident along the crystal axis — Cont. 

b, Crystals No. 5, 4, 3, and 1. 
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The statistical fluctuations of the Monte Carlo 
calculations, as evidenced by the irregularities of 
the histograms, are much worse for high than for 
low incident energies. It is noticeable that quite 
similar local irregularities occur for histograms 
pertaining to the same energy but different crystals. 
This phenomenon is due to the use of the same groups 
of photon histories for all the crystals. 

3.4. Other Source Geometries 

Monte Carlo calculations of the response function 
for any desired source geometry can easily be 
accomplished through the proper choice of the initial 
conditions, photon positions and directions, in the 
random sampling process. Sample results of such 
calculations for 0.661 Mev radiation are presented 
in figures 2, 3, and 4. 

a. Broad-Beam Radiation 

In figure 2 the ratio Pbb/p is plotted against r/R, 
where p BB is the photofraction for a broad cylindrical 
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Figure 2. Comparison of the photofractions for broad bea?ns, 
and for collimated radiation incident along the crystal axis. 
The ratio p BB /p is plotted as a function of r/R, where r is the radius of the broad 
beam, and R that of the crystal. The source energy is 0.661 Mev. 
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Figure 3. Comparison of the photofractions for point-isotropic 
sources, and for collimated radiation incident along the 
crystal axis. 

The ratio p IS /p is plotted as a function of tan-i (D/R), where D is the distance 
from the crystal to the point source (located on the crystal axis), and R is the 
crystal radius. 



beam with cross-sectional radius r, p is the photo- 
fraction for a collimated source, and R is the crystal 
radius. Three curves are shown, for crystals No. 3, 
large; No. 6, medium-sized; and No. 9, small. 

b. Point-isotropic Source 

In figure 3, the ratio pi 8 /p is plotted against 
tan -1 R/D where p IS is the photofraction for a point- 
isotropic source located at a distance D from the 
crystal on the cylinder axis. The curve pertains to 
crystal No. 6. The end-point of this curve, for 
ta,n~ 1 R/D=0°, was obtained from the corresponding 
broad-beam result. 

c. Off-Axis Collimated Beam 

In figure 4, the shape of the response function is 
shown for collimated pencils of radiation incident 
at various distances s from the crystal axis, for 
crystals No. 4 and No. 1. The normalization and 
other details are the same as in figure 1 ; the photo- 
fractions are also shown. 
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Figure 4. Response function of Nal(Tl) scintillation crystals 
for collimated radiation incident off center at a distance s from 
the crystal axis. 

The curves are normalized so as to correspond to one incident photon, the area 
under them being Y(E)[l—p(E)]. The curves for different values of s have been 
shifted apart through multiplication by indicated scale factors 10a. Only the 
"Compton tails" of the spectra are shown; the photopeaks are indicated by 
arrovvs, and the values of the photofractions are shown beside the corresponding 
curves. 
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4. Discussion 
4.1. Comparison With Other Calculations 

Our results for the value of the photofraction for 
collimated radiation incident on small crystals are 
in good agreement with those obtained by Maeder 
et al. For example, at an energy of 0.661 Mev 
one finds: 



p( Maeder) _ 
p (table 2). 


Crystal number 


9 


8 


7 


6 


5 


4 


0. 246 
. 243 


0. 368 
.377 


0. 418 
. 442 


0. 471 
.481 


0. 507 
.508 


0. 546 
. 542 



At 0.279 Mev and 1.17 Mev the agreement is 
almost as good (with the exception of a 15 percent 
discrepancy for the smallest crystal (No. 9) at 
1.17 Mev). 

No other results comparable to our photofractions 
for large crystals No. 1, 2, and 3 exist. Maeder et 
al. did, however, attempt to compute the photo- 
fraction for the limiting case of an infinitely large 
crystal (i. e., a semi-infinite medium), calculating 
the first two interactions of the incident photon 
exactly, and estimating the effect of the subsequent 
interactions. We have calculated the photofraction 
for this case by extrapolation based on the empirical 
formula (12) and arrive at the following comparison: 



Infinitely large crystal 


Energy (Mev) 


0. 279 


0. 661 


1. 17 


p( Monte Carlo) 


0. 975 

. 975 


0.940 

. 876 


0. 980 
. 802 


p( Maeder) 



The "orders of scattering" approach, in the limit 
of very large crystals, agrees with the more exact 
calculation only at 0.279 Mev, where multiple inter- 
actions are unlikely because of the large gamma-ray 
absorption cross section, but leads to a significant 
underestimate of the photofraction at higher energies. 

Our albedo values 1-p are in good agreement with 
the results of a calculation by Hay ward and Hubble. 13 
Note that the albedo is greater at 0.6G1 Mev than at 
0.279 or 1.17 Mev. There is also experimental evi- 
dence for the occurrence of this maximum at an 
intermediate energy. 14 

The pictures of the response function given by 
Maeder et al. are similar in appearance to ours, but 
not detailed enough for an accurate comparison in 
regard to shape. 

4.2. Comparison With Experimental Results 

For small crystals a detailed comparison with ex- 
periments has been made by Maeder et al., which 

13 E. Hayward and J. Hubbell, Phys. Rev. 93, 955 (1954). 

14 H. W. Koch (private communication). 



can equally well be applied to our results. It indi- 
cates fair agreement between theory and experiment; 
but there is a consistent tendency for the experi- 
mental photofractions to be somewhat lower than 
the theoretical values. This discrepancy may be 
reduced in part if one takes into account energy 
losses due to electron escape from the crystal. 15 

Measurements of Foote and Koch (see footnote 2) 
indicate for collimated radiation incident on a large 
crvstal (K—2.5 in., L=4 in.) photofractions p=0.S 
at"0.661 Mev, and 0.64 for cobalt-60 radiation, while 
the corresponding theoretical values are p= 0.821, and 
0.675 respectively. 16 

Further results for a large crystal (R=2 in., 
L=4 in.) have recently been reported by Kreger (see 
footnote 6). His incident radiation was a broad 
beam with cross-sectional radius r=0.25 R. Hence 
the comparison must be made for the broad-beam 
photofraction p BB . We have so far only calculated 
this correction at an energy of 0.661 Mev, but esti- 
mate that with an error of only a few percent tin' 
same fractional correction can also be applied at the 
other energies of Kreger's experiment. 



p (Kreger) _ . 

p (present results) _ 


Energy (Mev) 


0. 279 

0. 93 

(.90) 


0. 661 


1. 17 


1. 32 


2. 02 

0. 47 

(• 50) 


0. 725 

. 74 


0. 57 
(• 60) 


0. 54 

(.58) 



The agreement is seen to be quite good at 0.661 Mev, 
where the exact value of the correction is known, 
and fairly good elsewhere. 

Summing up the entire experimental evidence, 
both for small and large crystals, we can state the 
experimental photofraction is invariably below the 
theoretical values, by a smaller or greater amount. 
One is led to the conclusion that in the experiments 
there is always background radiation present (such 
as back-scattered radiation from the radiation source, 
and from the material surrounding the detector) 
which would be counted in the tail end of the pulse- 
height spectrum and would thus tend to depress the 
value of the photofraction. 



The authors thank Dr. Evans Hayward for several 
stimulating and enlightening discussions, and Miss 
Mary Orr for her assistance with the hand com- 
putations. 

5. Appendix: Escape of Secondary 
Radiation 

5.1. Effect on the Photofraction 

To estimate the effect of the escape of secondary 
radiation on the value of the photofraction p, we must 
determine the probability that an absorption is 

is K. Liden and N. Starfelt, Arkiv Fysik 7, 428 (1954). 

_ Y(1.17)p(1.17) + yq.33)p(1.33) 



is The photofraction for cobalt-60 radiation, p = - 



r(1.17)+F(1.33) 
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"spurious/ ' i. e., that it results in, or is preceded by, 
the escape of at least one secondary photon from 
the crystal (with an energy greater than 50 kev). 17 

We shall ignore the displacement of charged 
particles in the crystal. 18 Bremsstrahlung and 
annihilation quanta are thus assumed to be emitted 
at the point of absorption of the primary photon. 
Moreover, the emission is assumed to be isotropic. 
First, attention is concentrated on the spurious 
absorptions which occur in the first interaction of the 
incident photon (i. e., prior to any Compton scatter- 
ings), and return later to the relatively unimportant 
spurious absorptions resulting from other inter- 
actions. Secondary radiation is therefore considered 
as originating from an exponentially distributed line 
source on the axis of the crystal (on the assumption 
that the incident radiation is collimated). 

The following notation will be used for the various 
attenuation coefficients: mph = photoelectric absorp- 
tion coefficient, /z p =pair production absorption co- 
efficient, /ulc = Compton scattering attenuation coeffi- 
cient; /z=/iph-[-Mp+mc== total attenuation coefficient. 

Annihilation Radiation. If an absorption of a 
primary photon of energy E results in pair formation 
and subsequent annihilation of the positron, the 
probability that at least one of the two annihilation 
quanta will escape from a cylindrical crystal (with 
radius R and length L) is given, approximately, by 
the following expression: 



Vo= 



iW) 



"l-exp [~fx(E)L] 



s: 



dxe-" (E)x 



d cos d{l~{l-e- s ^) (l-e~ s ^)pfpf} (15) 

where /x / = J u(mc 2 ), and Si and s 2 are the distances 
traveled in the crystal by the two annihilation quanta 
(emitted in opposite directions). 

If Q<x<L/2 

cos 0i <cos < 1 Sx^x/cos s 2 =(L—x)/cos 

cos 6 2 < cos < cos 6i Si=x/cos s 2 =R/sin 

O<cos0<cos0 2 s l =R[sin0 s 2 =R/sm0 

If L/2<x<L 

cos 2 < cos 0<1 s 1 =xfcos s 2 =(L—x)/cos 

cos 0i < cos 0<cos 2 s^R/sin s 2 =(L—x)/cos 
O<cos0<cos0! s 1 =R/sm0 s 2 =R/sin0 

(16) 



where 
cos 0i= 



M2£p-*- ['+(I)T' 



17 The limit of 50 kev was selected because the same low-energy cutoff was also 
used in the main Monte Carlo calculation. 

i 8 The ranges for 0.3, 1.0, and 4.0 Mev electrons in Nal are approximately 0.02, 
0.17, and 0.861 cm, respectively. Hence only at the highest energy (4.45 Mev) 
and for the two smallest of the crystals considered (No. 8 and 9) could the neglect 
of electron displacement lead to a significant overestimate of the photofraction. 



p* and p 2 are photofractions for the two annihilation 
quanta, which — we believe — can be accurately ap- 
proximated by the photofraction p for an energy of 
0.511 Mev, as obtained by interpolation from table 2. 
Bremsstrahlung. The probability of escape of a 
bremsstrahlung quantum of energy E' can be ex- 
pressed by a formula similar to (15): 

^ l{fL) 1-exp [-»(E)L] Jo dXe 

\ f 1 dcos0 {e- s » {E ^ + [l-e- s ^ E,) ][l-p*(E)}} 

(17) 

where we again set p*(E)=j)(E), and the path 
length s is defined as follows: 



cos 2 <^cos 0<1; s= 



x 



cos 



- cos 0i <C cos < cos 2 ; s=— — - 

sm0 



— 1 <cos 0<— cos Oi] 



x — L 

cos 



(18) 



(cos 0i and cos 2 have the same meaning as above). 
Let <t>(Ej E')dE f denote the probability, per unit 
path length, for an electron or positron of energy E 
to produce a photon in the energy internal (E' f 
E'+dE'), 19 and let dE/ds be the average energy loss 
per unit path length of an electron or positron due 
to ionizing collisions. 20 The probability of an escape 
of a bremsstrahlung quantum produced by an 
electron or positron of energy E in the course of 
slowing down is 



P2(E)- 



~R(E) 



(E)Jz 



dE f 



J*'dE» ^^)/-^(^0] Pi(E"), (19) 

where R(E) is the electron or positron range, and 
E=5Q kev is the cutoff. 

Let \//(k)dk be the probability that one of the 
members of the pair has acquired initially a frac- 
tion between k and k-\-dk of the total available 
kinetic energy. 20 The probability of the escape of 
at least one photon in the form of bremsstahlung 
emitted by a pair of total energy E while slowing 
down is 

P,(E)= P dk+(k)[P 2 (kE)+P 2 (E-kE) 

-P 2 {kE\P 2 {E-kE)]. (20) 



i fl H. Bethe and W. Heitler, Proc. Roy. Soc. [A] 146, 83 (1934). 
20 H. Bethe, Handbuch der Physik, 24, 519 (1933). 



364 



The above formula assumes that each member of 
the pair can emit at most one photon .with energy 
greater than 50 kev in the course of slowing down. 
We believe that the effects of plum! photon produc- 
tion can be safely neglected in the energy range of 
interest here. 

Combined Effect. The probability that the ab- 
sorption of a primary photon of energy E in a first 
interaction will be spurious due to the combined 
effects of bremsstrahlung and annihilation radiation 
is 



P,(E)= 



Mph 



Mph+Mp 



P 2 (E) 



+ -^—[Pz(E-2mc 2 )+P -P3(E-2rnc 2 )Po]. 

(21) 

In evaluating tins formula we found that the neg- 
lect of bremsstrahlung (i. e., setting P 2 and P 3 =0) 
would lower P 4 by only a few percent (a maximum 
of ~ 11% for crystal No. 9 at 4.45 Mev). The prob- 
ability, per incident photon, of a spurious absorption 
in the first interaction is 

P,[E)=[Y(E)lix{E)] [^(E)+^(E)]P A (E). (22) 

Next we must determine a corresponding probability 
of spurious absorption in a second interaction. One 
contribution to tins probability is similar in origin to 
(21), and can be written: 



P*(E) = 



-[Y(E)/»(E)]j } 



ii(E f ) 



k(E,E')Q(E')Pf(E'), 



(23) 



where k(E,E f ) is the Klein-Nishina differential co- 
efficient for Compton scattering with energy change 
from E to E f ; Q is the probability that the scattered 
incident photon will not escape from the crystal with- 
out further interaction 21 ; and Pj(E) is the probabil- 
ity that a given absorption event is spurious, and is 
analogous to P 4 (E), but should take into account 
the fact that the source of secondary radiation is no 
longer concentrated on the crystal axis. We have 
ignored this circumstance, because the evaluation of 
an exact expression for P*(E) would be exceedingly 
complicated, and have set P%(E) = P 4 (E). It is 
hoped that this will introduce only a minor error, 
for the following reasons: For small crystals, such a 
geometrical error might be serious, but a second inter- 
action is quite rare; for large crystals, in which multi- 
ple interactions are more common, the geometrical 

2i An expression for Q has been given by Maeder et al. (footnote 3, eq 44, 45, 46;) 
their Ho is our YQ. 



error is much less significant. In any case, the 
contribution of P 6 to the change in the pnotofraction 
will be small, so that a rough-and-ready treatment is 
in order. 

There is also the possibility thai mi absorption 
taking place in a second interaction is spurious be- 
cause it was preceded by the escape of a brems- 
strahlung quantum from a. Compton recoil electron. 
The probability for this to happen is 



P 7 (E)=[Y(E)/n(E)]f*dE 



, »ME')+ n*(E') 



n(E') 



k(E 1 E')Q(E')P 1 (E-E'). 

Thus a spurious interaction will occur 
probability 



(24) 
with 



P 8 (E)=P 6 (E)+P 7 (E)~P Q (E)P 7 (E) (25) 

per incident photon in a second interaction; and with 
probability 

P t (E)=P s (E)+P a (E)-P 5 {E)P a (E) (26) 

in a first or second interaction. At the highest 
energy and for the largest crystal considered, the 
value of P s was only 15 percent of P 5 , and generally 
it was much smaller. Subsequent interactions will 
lower the photon energy further and make spurious 
absorptions even more unlikely. Hence we felt jus- 
tified in ending our calculations with the effects of 
the second interaction. The change in the photo- 
fraction is thus 



Ap=-P 9 (E)Y(E). 



(27) 



5.2. Effect on the Shape of the Response Function 

In considering the detailed shape of the response 
function we have calculated only the effects of an- 
nihilation radiation from a first interaction. Other 
escaping secondary radiation was taken into account 
only insofar as it shifted pulses from the photopeak 
into tail of the response function (i. e., by an appro- 
priate renormalization of K(E,E')). 

None, or one, or both of the annihilation quanta 
resulting from a pair production can escape from the 
crystal without further interaction. The probabil- 
ities for these contingencies can be obtained by 
replacing the integrand in the inner integral in (15) 
by the expressions 



(a) (no escape) : (£— e~W)Q— g-V) 

(b) (one escape) : e~ s ^'(l — e~ s ^)-\-e~ s ^' 

(l—e-W) 

(c) (two escapes): e~ s ^'e~ s ^'. 



(28) 
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In case (a) the energy E' left in the crystal is the 
sum of three terms; two of them are each distributed 
as 

L(mc 2 )/Y(E)=K( mc 2 E r ) +p(mc 2 )5(E' - mc 2 ) (29) 

the third contribution is a constant amount of 
energy E— 2mc 2 . 

In case (b), E f is the sum of two terms, one dis- 
tributed according to (29), the other a constant 



equal to E— 2mc 2 . In case (c) £ / = £-2mc 2 . The 
correction of the response function was performed 
by shifting the appropriate number of pulses from 
the photopeak to the "tail" of the pulse height 
spectrum, distributing them according to the dis- 
tribution functions described above weighted ac- 
cording to (28) . 



Washington, November 4, 1955. 
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